Diquark correlations in baryon spectroscopy and holographic QCD 
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We introduce an improved mass formula for the nucleon and delta resonances and show how it 
emerges from AdS/QCD in a straightforward extension of the "metric soft wall" gravity dual. The 
resulting spectrum depends on just one adjustable parameter, characterizing confinement-induced 
IR deformations of the anti-de Sitter (AdS) metric, and on the fraction of "good" (i.e. maximally 
attractive) diquarks in the baryon's quark model wave function. Despite its simplicity, the predicted 
spectrum describes the masses of all 48 observed light-quark baryon states and their linear trajectory 
structure with unprecedented accuracy. 
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The arguably most prominent and pervasive struc- 
ture in the known hadron spectrum consists of (approxi- 
mately) linear Regge-type trajectories of equal slopes on 
which the square masses M 2 of excited states organize 
themselves in Chew-Frautschi plots, i.e. as a function of 
either the angular momentum L (or total spin J) or the 
radial excitation level N. The QCD-based understand- 
ing of these trajectories and their relation to linear quark 
confinement remains one of the pre-eminent challenges of 
strong-interaction physics. 

In this Letter we propose a mass formula for the trajec- 
tories in the light-quark baryon (i.e. nucleon and delta) 
sector which improves upon the recent AdS/QCD |l| 
"metric soft-wall" (ms) prediction Q 



N ■ 



(1) 



and show how it emerges from a straightforward exten- 
sion of the ms gravity dual. More specifically, while Eq. 
([1]) works very well in the A sector (i.e. all observed 
A* resonance states lie within errors on the predicted 
trajectory Q whose empirical slope is used to determine 
the AdSs deformation parameter A = 0.52 GeV), we will 
show that its predictions for the nucleon excitation spec- 
trum can be substantially improved by adding the cor- 
rection 



AM L = - 2 ( M A ~ M n) ^ 



(2) 



which solely depends on the resonances' diquark con- 
tent. The latter enters through the good diquark frac- 
tion K g d which we define as the fraction of the space-spin- 
flavor baryon wavefunction in the SU(2) S <g> SU(3)^ prod- 
uct representation which contains the most attractive, 
i.e. the "good" (0 + color-antitriplet iso-singlet) diquark. 
This implies in particular K gc j = for all A and spin-3/2 
nucleon resonances, K g d = l/4 for spin- 1/2 nucleons in 
the 70su(6) representation and K g d = 1/2 for nucleons in 
the 56su(6) multiplet. The mass correction © renders 
the long-suggested 0| and recently reemphasized H 0] 
importance of attractive diquark correlations in baryon 



structure explicit. Interest in such correlations currently 
extends to phenomena ranging from exotic hadrons 6 1 to 
color-superconducting phases of ultradense matter [7 1 . 

In order to compare Eqs. (jTJ) , to experimental data, 
one needs to assign intrinsic orbital and spin angular mo- 
menta L and S to the observed states. At first glance this 
seems unfeasible since only total spin J and parity P are 
measured. However, the data show regularities which 
bring additional order into the spectrum (see Table []}. 
The three-body dynamics of the quark model implies that 
the number of expected states grows dramatically with 
increasing excitation energy. We will discuss these states 
in the non-relativistic harmonic oscillator (h.o.) approxi- 
mation. Although the h.o. states may mix, we will show 
that the leading L and S configurations can nevertheless 
be determined. In terms of the Jacobi coordinates p and 
A, the two h.os. exhibit orbital (l„, l\) and radial (n„, n\) 



excitations with L= l +l\ and N= 



-n\, and the en- 



ergy levels develop a shell structure which manifests itself 
experimentally in a series of resonance regions. 

The Particle Data Group lists 22 nucleon and 22 A 
resonances with known J and P 8] . We use all of them 
except for the one-star A(1750) - which is special and 
requires a devoted discussion - and except for the two- 
star A(2000) with J p =5/2 + for which three mass values 
are quoted, 1724, 1752, and 2200 MeV. The long debated 
question whether the Roper resonance AT(1440) is the 
first radial excitation of the nucleon has recently found 
an affirmative answer in electro-production experiments 
9]. A(1600) with J p =3/2 + is the analogue state in the 
A spectrum; we asign L, N=0,1 to both. We further 
conjecture (A) that radial and orbital excitations cost 
the same energy, in line with the N + L dependence of 
Eq. (HJ). 

The two states iV(1535) and iV(152Q) with J p =l/2~, 
3/2-, and the three states AT(1650), N (1700), N (1675) 
with J p =l/2~, 3/2 _ , 5/2 ~ naturally form a spin dou- 
blet with L=l,5=l/2 and a £=1,5=3/2 triplet. Sym- 
metry arguments forbid a symmetric spin wave function 
for the lowest-mass A states with L=l, and indeed only 
the L=l, 5=1/2 spin doublet A(1620), A(1700) is ob- 
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TABLE I: Masses 
together with the 



and quantum numbers of N and A states, 
predictions of Eqs. (fTj) and (fl?jl . 



not have L = 5 (apart from a small admixture), rather 
it must have L = 3,5 = 3/2. The three other members 
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served. £=1,5=3/2 would form a triplet; symmetry re- 
quires excitation of the radial wave function to N=l. 
These states are seen as A(1900), A(1940), and A(1930) 
having J p =l/2~, 3/2~, 5/2~ quantum numbers. Their 
low mass - compared to quark model calculations - arises 
from conjecture (A). Furthermore, we conjecture (B) 
that all negative parity A excitations with 5=3/2 have 
N=l. 

In addition, there are three positive-parity nucleon 
resonances 7V(1900), 7V(1990), iV(2000) with quantum 
numbers J p = 3/2+, 5/2+, 7/2+, four A states A(1910), 
A(1920),A(1905),A(1950) with J p = 1/2+, 3/2+, 5/2+, 
7/2+, and no other close-by state with these quantum 
numbers. It is therefore natural to assign L=2 and spin 
5=3/2, yielding a A and a nucleon quartet; the latter 
quartet is completed by the recently proposed 7V(1880) 
with J p =l/2+ Q. 

The J p = 1/2+ 7V(2100) has no spin-parity partner; we 
interpret it as the third excitation in the series iV(940), 
7V(1440),iV(1710),iV(2100). The two states iV(2070), 
7V(2190) could form a L=3, 5=1/2 spin doublet, and 
the three states iV(2200), iV(2190), iV(2250) an incom- 
plete quartet where the lowest angular momentum state 
with J p = 3/2" is missing [lj|. The iV(2190) assignment 
is ambiguous. The J p =l/2~ 7V(2080) mus t have L=l 
and could be the second radial excitation of iV(1535). 
Its spin doublet partner is then ^(2090) as the second 
7V(1520) excitation, its isospin partner is A(2150). The 
isospin partner of iV(2190) is missing. First excitations 
have been suggested: 7V(1905) with J p =l/2" [ll|, [H 
and 7V(1960) with J p =3/2" [H, H Q- Their pre- 
dicted mass is 1.82 GeV. 

Due to its low mass, the J p = 9/2~ state 7V(2250) can- 
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with lower J are missing. A(2223) (suggested in 
and A(2200) with J p =5/2",7/2" are natural candi- 
dates for the L=3, 5=1/2 doublet. At a bit higher 
mass, A(2400) with J p =9/2~ is observed. According to 
our conjecture (B), it has L=3, 5=3/2, N=l. A(2350) 
with J p =5/2~ could be its spin-multiplet partner while 
J p = 3/2~, 7/2~ are missing. Conjecture (B) also entails 
the L = 5, 5 = 3/2, A^ = l assignment for the high-mass 
J p =13/2- resonance A(2750). At 2.4GeV, a nearly 
complete spin multiplet with L = 4, 5 = 3/2 is known: 
A(2390)(7/2+), A(2300)(9/2+), A(2420)(ll/2+). The 
assignments of 7V(2700) and A(2950), finally, are unam- 
biguous. 

Table[T]reveals the extent of the agreement between the 
48 measured mass values and Eqs. ([T]), ([2]). In fact, the 
latter considerably improve upon any dynamical quark 
model prediction of the full mass spectrum. This sup- 
ports the validity of both Eqs. |T]), ((2]) and of our quan- 
tum number assignments, and it raises the hope that not 
only Eq. |T]) but also the universal correction term @ 
may have a transparent origin in holographic QCD. In 
the following we will show that this is indeed the case. 
Since the metric soft wall flj is so far the only AdS/QCD 
dual which predicts linear square-mass trajectories for 
baryons (cf. Eq. (Q])), we take this background as our 
starting point. The IR deformed AdS5 line element 

ds 2 = gMN<lx M dx N — a 2 (z) {rj^ v dx^dx v — dz 2 ) (3) 

(where 77 is the four-dimensional Minkowski metric and 
z the conformal coordinate of the fifth dimension) of the 
metric soft wall in the baryon sector is characterized by 
the warp factor 



(4) 



(where R is the curvature radius, m[ ' the Dirac bulk 
mode mass and A the scale which governs the IR defor- 
mation). In this geometry the Dirac equation for the 
string modes (x, z) dual to nuclcons reads 




(d M + 2a( ms )- 1 cW ms) ) 



* (x,z) = 
. (5) 

(where ef.j = a5f f are fiinfbeins and T A Dirac matrices). 
By iteration, Eq. ([5]) can be cast into the Sturm-Liouville 
form 

_Q2 + yirns) ( J ^ = M (ms)2^(m S ) ^ (g) 

where the modes ip± (z) = (Aze - " 4 ) " <p± (z) describe 
the rescaled chiral components of the bulk spinor 

d 4 k 



* {x, z) 



— ikx 



\pf (z)P + 



+ tp W (z) P_ $W (fc) (7) 



,(4) 
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(P± = (l ± 7 5 ) /2) and the boundary spinor solves 
the d = 4 Dirac equation (j a k a - \k\) *W (k) = 0. The 
metric soft-wall potential is [2J 



V 



(ms) 
,C°> - 



(z)=±d z {a^m^)+a^ 2 mf )2 (8) 



(ms) 



can be found 



analytically Q. After implementing the boundary be- 

havior <ys + (z) — > z T required by the AdS/CFT dictio- 
nary for baryon interpolators of twist dimension f , which 
amounts to replacing m[ ' R — s-f — 2 = L + 1 [l6| , the 
eigenvalue spectrum (fTJ) emerges. 

As mentioned above, this spectrum agrees very well 
with data in the A sector but can be substantially im- 
proved in the nucleon sector by adding the correction 
term ([2]) which depends on the diquark content of the 
baryon resonances. At first it might seem impossible 
to gain holographic access to this diquark information 
because diquarks and their operators are gauge depen- 
dent whereas only gauge-invariant operators have well- 
defined duals. However, diquark information may enter 
indirectly through gauge-invariant baryon interpolating 
fields whose general form [ItJ 



r\t (x) = 2 [r) pd (x) + tr) sd (x)} 



(9) 



at leading twist (i.e. with the minimal scaling dimen- 
sion 9/2) contains a pseudoscalar diquark in i] p d = 
£abc (u^Cdb) j 5 u c and a "good" scalar diquark in ?7 s d = 
£ahc (v^C-l^db) u c . (Higher-twist interpolators contain 
additional covariant derivatives.) 

The interpolators © are expected to have enhanced 
overlap with nucleon states of the corresponding diquark 
content and can thus be related to the n g d of these states 
[23l | . Since t — * 00 corresponds to maximal "good" di- 
quark content and t = — 1 (the "Ioffe current") to an 
exclusively "bad" diquark content, a simple approxima- 
tion to this relation may e.g. be 



t (« gd ) 



1 



2^gd — 1 



(10) 



The different diquark content of the baryon interpola- 
tors rjt can affect properties of their dual modes in sev- 
eral ways. For once, the AdS/CFT dictionary relates the 
mass 777,5 of a dual string mode to the scaling or twist 
dimension of the corresponding interpolator. While all 
equal-twist baryon interpolators have the same classi- 
cal twist dimension f = L + 3 and hence correspond 
to dual modes of equal mass m§' = (L + 1) /R [l6j], 
their anomalous dimensions will generally depend on 
t and therefore induce K g d dependent mass corrections 
A7715 (ftgd) ■ Such corrections could also result from K g d 
dependent couplings of the Dirac modes to condensate- 
related spin-0 bulk fields. Since at present no reliable 
QCD information on nonperturbative anomalous dimen- 
sions exist, we will treat A7775 as an adjustable parameter 



to be determined in bottom-up fashion from the observed 
baryon spectrum (as similarly suggested in Ref. [2lj). 

We are now going to show that such a bulk mass cor- 
rection A7775 (K g d) indeed results in the desired spectral 
correction ((2|) and minimally just requires a small addi- 
tional IR deformation of the soft wall metric ([3]). To this 
end, we first observe that a shifted dual mode mass 



(0) 1 A 
7775 = 7775 + A77I5 = 



L + Am 5 i? + 1 
R 



with the corresponding IR deformation 



. . R/ A 2 z 2 

a (z) = — H 7 

w z \ L + Am 5 R + 1 



(11) 



(12) 



of the warp factor Q turns the Sturm-Liouville potential 
© into 

V L , ± (z) = V^ +Am5 ± (z) = Vgt^ (*) • (13) 

Hence the eigenfunctions ip± in this potential can be ob- 
tained from the analytical solutions ip 1 ^ 1 ^ ± of Ref. 

as ipN,L, Kgd ,± (z) = i>N$ + Am B R,± ( z ) ( wnich are wel1 " 
defined for Am$R > —3/2). Finally, after specifying 



Am 5 (Kgd) 



AMI 
4X 2 R 



(14) 



(so that Am^R > —0.7 and 777,5 > 0) the eigenvalue spec- 
trum ([1]) indeed turns into 



M 2 N , L 



M 



(ms)2 

N,L+Am 5 R 



4A \ N + L 



am: 



The above derivation also sheds new light on the physics 
behind this spectrum. Indeed, it shows that the dual 
modes ipN,L,K sd ,± corresponding to larger K g d (with iden- 
tical N, L) begin to feel the soft wall at smaller z (cf. Eq. 
(|12p) and therefore extend less into the fifth dimension. 
This reflects the additional attraction and translates into 
a smaller size of baryon states with larger K g a- 

A more elaborate alternative for generating the spec- 
trum (fT"5l) would take the RG flow of the anomalous di- 



mensions (or additional bulk fields) into account [24 1. 
According to the gauge/gravity correspondence, this 
renormalization-scale dependence translates into a z de- 
pendent Amc,(z). Since the above adaptation of the 
eigenfunctions then ceases to work, one may instead at- 
tempt to generate the mass correction directly in the 
mode potential, 



V L , ± (z) = (z) 



am: 



(16) 



while both the warp factor (01 and the eigenfunctions 
of the metric soft wall remain unchanged. Since the z 



4 



dependent mass correction Ato5.± (z) = m± (z) results 
in the modification 



AV± (z) =±(a 



(ms) m± 



a 



(ms)2 ~ 



m± ( 2m 



(0) 



m± 



(17) 

(the prime denotes d z ) of the soft-wall mode potential 
©, a 2 (and L) independent shift AV± = AM^ gd is 
obtained when the m± (z) solve the nonlinear differential 
equations 







1 , s fl( ms )' 



m± 



(ms) 



(18) 



General solutions of these equations indeed exist and, 



subject to the boundary condition rh± (z) — ► 0, read 



m± (z) = 



4 (L + 1) Am 5 X 2 z : 
L + 1 + X 2 z 2 



f±(z) 



(19) 



with 



/+(*) 



1 1 F 1 (1-Am 5 i?,i+|,-A 2 z 2 ) 



2L + 3 1F1 (-Am 5 fl, L + 
1 iFi(l + Am 5 fl,-Z, 
2L+1 1F1 (Am 5 i?, -L - 



U 2 ^ 



,A 2 ; 



-,(21) 



where the constant Ams is given by Eq. (fT4")l and 
1F1 (a, 6, x) are Kummer functions [22| . Inserted into Eq. 
(fl~7|) . the solutions (fl"9]) reproduce Eq. (fl~6)) and hence the 
spectrum (|15p . However, it is not immediately obvious 
how the chirality dependence of the mass functions (fT^)) 
can arise directly from a modified Dirac mass term. 

To summarize, we have presented a light-quark baryon 
mass formula with only one adjustable parameter (re- 
lated to the string tension of linear quark confinement) 
which reproduces the masses of all 48 observed nucleon 
and A resonances with far better accuracy than e.g. 
quark model predictions (although the latter depend on 
a substantially larger number of parameters). In addi- 
tion, our spectrum relates the trajectory slopes to the A 
ground-state mass and reveals a strikingly systematic role 
of the good diquark fraction n g d in baryon spectroscopy. 

We have furthermore shown how this spectrum 
emerges in holographic QCD, namely by means of string- 
mode mass corrections in the metric soft wall background 
which can be naturally traced to the varying diquark con- 
tent of the QCD baryon interpolators. The latter reflects 
itself e.g. in different anomalous dimensions and poten- 
tially in diquark-contcnt dependent couplings of the dual 
modes to condensate- related spin-0 bulk fields. (It would 
be interesting to gain independent quantitative insight 
into these effects, e.g. by lattice calculations of the QCD 
baryon interpolators' anomalous dimensions in the in- 
frared.) Finally, our results indicate that baryon sizes 
decrease with increasing good-diquark content. 
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